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The Middle East respiratory syndrome coronavirus (MERS-CoV) is a newly identified pathogen able of 
human transmission that causes a mortality of almost 40%. As in the case of SARS-CoV, MERS virus 
lacking E protein represents a potential vaccine. In both cases, abolishment of channel activity may be 
a contributor to the attenuation observed in E-deleted viruses. Herein, we report that purified MERS- 
CoV E protein, like SARS-CoV E protein, is almost fully a-helical, has a single a-helical transmembrane 
domain, and forms pentameric ion channels in lipid bilayers. Based on these similarities, and the proposed 
involvement of channel activity as virulence factor in SARS-CoV E protein, MERS-CoV E protein may 
constitute a potential drug target. 


© 2015 Elsevier B.V. All rights reserved. 


1. Introduction 


In September 2012 the first case of a novel coronavirus (CoV) 
infection was reported (Zaki et al., 2012) and the etiological 
agent was named Middle East respiratory syndrome coronavi- 
rus (MERS-CoV) (de Groot et al., 2013). As of 23 May 2014, 635 
cases and 193 deaths have been confirmed (http://www.who.int/ 
csr/don/2014_05_23_mers/en/). An investigation on 47 cases 
between September 2012 and June 2013 revealed that MERS 
patients display similar clinical features to severe acute respira- 
tory syndrome (SARS) virus infection, but progress more rapidly 
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towards respiratory failure (Assiri et al., 2013). In addition, a major- 
ity of MERS patients (estimated at 75%) has at least one co-morbid 
condition, in contrast to only 10-30% in SARS patients (Assiri et al., 
2013; Drosten, 2013; Raj et al., 2014). Mortality rate in MERS-CoV- 
infected individuals is about 40%, compared to 10% during the SARS 
outbreak (http://www.who.int/csr/sarsarchive/2003_05_07a/en/). 
Human-to-human transmission has been reported, although it is 
still limited and inefficient (Raj et al., 2014). Although there are 
still no effective treatments, a mouse model suitable for studying 
MERS-CoV infection has been engineered (Zhao et al., 2014). The 
genomic sequence of MERS-CoV is closely related to bat CoVs HKU4 
and HKUS5, which belong to clade c of B-coronaviruses (Annan et al., 
2013; Drexler et al., 2014). Like SARS-CoV, MERS-CoV is thought to 
have crossed over from bats to humans, although MERS-CoV used 
camels as an intermediate species (Haagmans et al., 2014). 

The structural proteins in coronaviruses are N (nucleocapsid) 
and integral membrane proteins S (spike), M (membrane) and E 
(envelope). The S protein is involved in host recognition and entry 
into the cell. Recently, the crystal structure of the fusion core (Gao 
et al., 2013) and the receptor binding domain (Chen et al., 2013) of 
MERS S protein have been solved. The M protein is the most abun- 
dant protein in the lipid envelope of virions and defines their shape 
(Barcena et al., 2009; Neuman et al., 2011), although high resolu- 
tion structural data is still not available. The E protein is scarcely 
present in the virion, but abundantly expressed in the infected 
cell, mainly distributed in intracellular membranes between ER and 
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Golgi compartments (Lim and Liu, 2001; Nalet al.,2005; Raamsman 
et al., 2000). In those sites, it participates in virus assembly, budding 
and intracellular trafficking. In the case of SARS-CoV, the E protein 
accumulates at the ERGIC (Nieto-Torres et al., 2011). 

CoV E proteins are 76-109 amino acids long and predicted to 
have at least one a-helical transmembrane (TM) domain based 
on prediction algorithms (Surya et al., 2013). The extramem- 
brane C-terminal domain has some hydrophobic character, and 
could represent in some cases a second TM domain, e.g., in infec- 
tious bronchitis virus (IBV) and in murine hepatitis virus (MHV) 
E proteins, although this has not been confirmed experimentally. 
SARS-CoV E protein has been shown unequivocally to have just 
one TM domain (Li et al., 2014), and is the only E protein for 
which structural data is available (Li et al., 2014; Pervushin et al., 
2009). E proteins have been shown to form ion channels in syn- 
thetic membranes (Torres et al., 2007; Wilson et al., 2006), with 
reported conductance of 20-30 pS. However, characterization with 
properly purified full length and transmembrane peptides of SARS- 
CoV E protein has produced much higher (>10 times) conductances 
(Baguena et al., 2012; Verdia-Baguena et al., 2012, 2013a,b). In 
SARS-CoV, this channel activity has been linked to activation of 
caspase-1 and maturation of pro-inflammatory IL-1 during infec- 
tion (Nieto-Torres et al., 2014), although a detailed mechanism to 
explain this functional link is still not available. 

Removal of E protein is in general deleterious to coronaviruses, 
although the severity of the damage is strongly virus-dependent. 
For example, deletion of E gene in transmissible gastroenteritis 
virus (TGEV) results in blockage of virus trafficking in the cellu- 
lar secretory pathway and prevention of virus maturation (Ortego 
et al., 2002). Mutations at the C-terminal region of MHV E protein 
results in impaired virion assembly and maturation (Fischer et al., 
1998). SARS-CoV lacking the E gene is attenuated (DeDiego et al., 
2007), and MERS-CovV lacking the E gene cannot propagate, despite 
showing similar viral RNA level to wild-type virus (Almazan et al., 
2013). Based on the attenuating effect of E protein removal, a SARS- 
CoV lacking E gene is currently being studied as potential vaccine 
(Enjuanes et al., 2011; Fett et al., 2013; Netland et al., 2010), anda 
similar E-deletion approach has been suggested for controlling the 
MERS outbreak to obtain a virus that is replication-competent, but 
propagation-defective (Almazan et al., 2013). 

The only coronavirus E protein structurally characterized in 
some detail is the one in SARS-CoV, as a transmembrane domain 
(Pervushin et al., 2009; Torres et al., 2006) and latter as a truncated 
form (residues 8-65), obtained using a 20 residue N-terminal tag 
that included a poly-His tail (Li et al., 2014). Herein, we report for 
the first time the expression, purification, and preliminary struc- 
tural study on full-length 82 residue MERS-CovV E protein without 
any tag. This purified protein allowed its biophysical characteri- 
zation to conclude that it has striking similarities to SARS-CoV E 
protein. 


2. Materials and methods 
2.1. Protein expression and purification 


The DNA sequence encoding E protein from MERS-CoV-EMC12 
(GenBank accession number JX869059) was fused N-terminally 
to a 6-His-MBP tag, followed by a TEV protease cleavage site, to 
form the construct 6-His-MBP-TEV-E. This was subcloned into pTB- 
MalE for expression. To investigate the effect of the four native 
cysteines (Fig. 1A), in addition to wild type MERS-CoV E, three addi- 
tional Cys-to-Ala mutants were generated, at positions 23 and 30 
(Eaacc), 40 and 43 (Eccaa), and 23, 30, 40 and 43 (Eaaaa). Expression 
and purification of non labeled and !°N-labeled Eqaaa Was per- 
formed following the protocol described previously for MBP-tagged 
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Fig. 1. Purification and secondary structure of MERS-CoV E protein. (A) Amino acid 
sequence of full-length E protein used in this study with predicted TM region in bold 
and underlined. The three N-terminal amino acids, SNA (gray), result from MBP tag 
removal by TEV protease cleavage. Positions of Cys residues, which were mutated 
into alanine, are indicated by ‘A’; (B) SDS-NuPAGE of purified Cys-less mutant, Eaqaa. 
The protein band migrated approximately at its calculated molecular weight of 
9.6 kDa; (C) infrared amide I region of MERS-CoV E in hydrated POPC lipid bilay- 
ers (solid line) and its corresponding Fourier self-deconvolved spectrum (dashed 
line). 


full-length SARS-COV E protein (Li et al., 2014), except that the 
elution buffer contained 1.6 mM n-dodecyl B-p-maltoside (DDM, 
Calbiochem). 


2.2. Gel electrophoresis 


For SDS-NuPAGE, samples were run in 4-12% NuPAGE® Bis-Tris 
gels (Invitrogen) with NuPAGE® MES SDS running buffer and 
stained with SimplyBlue™ SafeStain (Invitrogen) according to the 
manufacturer’s protocol. Blue-native PAGE (BN-PAGE) was per- 
formed as described previously (Gan et al., 2011). Lyophilized 
peptide was solubilized to 0.1 mM in a sample buffer containing 
50 mM DHPC (Avanti Polar Lipids), 5, 25 and 50 mM DPC (Avanti 
Polar Lipids), or 21 mM LMPG (Anatrace). Escherichia coli aquaporin 
Z(AqpZ) in20 mM SDS was included as additional molecular weight 
marker. 


2.3. Fourier-transform infrared spectroscopy 


Sample preparation, ATR data collection and H/D exchange were 
performed as described previously (Torres et al., 2006) on a Nico- 
let Nexus spectrometer (Madison, USA). Briefly, the peptide was 
incorporated in multilamellar liposomes by dissolving a dry mix- 
ture of peptide and POPC (Avanti Polar Lipids) or DMPC (Avanti 
Polar Lipids) at a 1:50 molar ratio. Fourier self-deconvolution was 
performed on the amide | region using FWHH of 20cm7! and a 
narrowing factor k of 1.5 (Kauppinen et al., 1981). 


2.4, Analytical ultracentrifugation 


Sedimentation equilibrium data were collected at 20°C ina 
Beckman XL-I analytical ultracentrifuge monitoring the absorbance 
at 280nm. Lyophilized peptide was dissolved to OD289 of 0.3, 
0.5, and 0.8 (12mm pathlength) in 50mM Tris buffer (pH 7.3), 
100mM NaCl, 5mM C14-betaine and 32.3% D20. The samples 
were centrifuged in six-channel charcoal-filled Epon centerpieces 
with quartz windows. Radial distribution profiles were acquired at 
18,500, 22,700, and 27,800 rpm after reaching equilibrium. The lat- 
ter was tested by HeteroAnalysis. Data were processed and fitted to 
several monomer-N-mer models in SEDFIT and SEDPHAT (Schuck, 
2003). 
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Fig. 2. HSQC spectrum of !5N-labeled MERS-CoV E protein. (A) ['H-!5N]-HSQC spectrum of 0.2 mM MERS-CoV E in 200 mM SDS micelles dissolved in H20; (B) same as (A), 
after exposure to 99% D20; (C) alignment of E proteins from MERS-CoV and MHV. The position of the four Cys residues in MERS-CoV E is indicated. 


2.5. Conductance measurements 


The experiments reported here have been carried out 
in single ion channels reconstituted on planar membranes. 
Planar bilayer membranes were formed from two mono- 
layers prepared from a lipid mixture containing 1,2-dioleoyl- 
sn-glycero-3-phpsphocoline (DOPC), 1,2-di-(9Z-octadecenoyl)-sn- 
glycero-3-phospho-.-serine (DOPS) and 1,2-dioleoyl-sn-glycero-3- 
phosphoethanolamine (DOPE) at a 3:1:1 (DOPC:DOPS:DOPE) in 
pentane on 70-90 pm diameter orifices in the 15 m-thick Teflon 
partition that separated two identical chambers (Bezrukov and 
Vodyanoy, 1993; Montal and Mueller, 1972). The orifices were 
pretreated with a 1% solution of hexadecane in pentane. Aque- 
ous solutions of KCl and NaCl were buffered with 5mM HEPES 
at pH 6. All measurements were performed at room tempera- 
ture. Single channel insertion was achieved by adding 2-3 pL of 
a 150 g/mL solution of protein in the buffer that contains acetoni- 
trile:isopropanol (40:60) at the cis side of the chamber. 

Anelectric potential was applied using Ag/AgCl electrodes in 2 M 
KCl, 1.5% agarose bridges assembled within standard 250 ml pipette 
tips. The potential was defined as positive when it was higher on 
the side of the protein addition (the cis side of the membrane cham- 
ber), whereas the trans side was set to ground. An Axopatch 200B 
amplifier (Molecular Devices, Sunnyvale, CA) in the voltage-clamp 
mode was used to measure the current and applied potential. The 
chamber and the head stage were isolated from external noise 
sources with a double metal screen (Amuneal Manufacturing Corp.., 
Philadelphia, PA). The single-channel conductance was obtained 
from the current measurement at applied potential of +100 mV in 
symmetrical salt solutions. 


2.6. NMR spectroscopy 


NMR experiments were performed at 318 K using an Avance-Il 
600 NMR spectrometer with a cryogenic probe. Sodium 2,2- 
dimethyl-2-silapentane-5-sulfonate (DSS) was used as the internal 


reference for 'H nuclei. The chemical shifts of !°N nuclei were cal- 
culated from the 'H chemical shifts. The NMR data were processed 
using TopSpin 3.1 (www.bruker-biospin.com) and analyzed using 
CARA (www.nmr.ch). To identify membrane-embedded residues, 
the NMR sample was lyophilized overnight and reconstituted in 
99% D20. Immediately after reconstitution, a 2D ['H-!°N]-HSQC 
spectrum was collected. 


3. Results and discussion 
3.1. Expression and purification of MERS-CoV E protein 


Full length MERS-CoV E protein was successfully expressed as 
a fusion protein with maltose binding protein (MBP) attached to 
the N-terminus of E. MBP was subsequently removed by TEV pro- 
tease, leaving behind additional 3 residues, SNA, at the N terminus 
of MERS-CoV E protein (Fig. 1A). Four mutants were tested (see Sec- 
tion 2), but only a Cys-null mutant (Eaaqaq ) could be obtained with 
sufficient purity, which was further characterized. This mutant is 
referred to hereafter as ‘MERS-CoV E protein’. The purity of this 
mutant was at least 90% after reverse-phase HPLC purification, as 
shown by gel SDS electrophoresis (Fig. 1B). When reconstituted in 
POPC lipid bilayers, the polypeptide was predominantly a-helical, 
as indicated by an infrared amide | band centered at 1655cm7! 
(Fig. 1C). These spectral features are also common to SARS-CoV E 
protein (Li et al., 2014), which suggests a similar overall secondary 
structure. 


3.2. MERS-CoV E has a single TM domain 


Screening of MERS-CoV E protein in micellar and bicellar envi- 
ronments identified SDS as the best detergent to achieve good peak 
dispersion in the two dimensions of an [!H-!5N]-HSQC spectrum 
(Fig. 2A), implying foldedness and lack of aggregation. We note 
that although this SDS environment will not allow the formation of 
oligomers (Fig. 1B), it is nevertheless likely to produce near-native 
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Fig. 3. Oligomeric state of MERS-CoV E protein. (A) Blue-native PAGE of 0.1 mM MERS-CoV E protein at the indicated concentrations of DPC, DHPC, and LMPG detergents. 
Oligomeric sizes are indicated by asterisks; (B) radial distribution profiles of the protein solubilized in C14-betaine micelles (open dots) and fitting to a monomer-pentamer 
self-association model (solid line) at the angular speeds shown. Residuals of the fitting are shown in the lower panel; (C) Chi-squared values corresponding to fitting to 
different monomer: n-mer models, from trimers to heptamers. The most optimal fit (monomer-pentamer) is indicated by an arrow; (D) homo-pentameric model of the TM 
domain of SARS-CoV E protein. The equivalent positions of the two Cys residues in MERS-CoV E protein, mutated in the present work, are highlighted (yellow, arrows) and 
shown at an interfacial orientation. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


secondary structure. This assumption is based on results obtained 
previously for SARS-CoV E protein (Li et al., 2014), where titration of 
the polypeptide SDS with increasing DPC, up toa 1:4SDS/DPC molar 
ratio, did not lead to changes in secondary structure, as monitored 
by !3C, chemical shifts. 

To determine if MERS-CoV E protein has one or two TM 
domains, the membrane-embedded residues were quantified using 
['H-!°N]-HSQC spectra in an HN-D20 exchange experiment. In 
99% D20, ~20 residues were found to be protected from exchange 
(Fig. 2B), which suggests a single a-helical TM domain. 

Early studies of N- or C-terminally tagged MHV E (Maeda 
et al., 2001) and SARS-CoV E (Yuan et al., 2006) proposed that 
both N and C-terminal ends were facing the cytoplasm, with the 
implication of two TM domains. However, purified full length SARS- 
CoV E protein has only one TM domain, as shown in a HN-D20 
exchange experiment (Li et al., 2014). In addition, the topology 
of untagged SARS-CoV E protein, both in SARS-infected cells and 
transfected cells shows that only one TM domain is possible, with 
C-terminal and N-terminal ends located cytoplasmically and lume- 
nally, respectively (Nieto-Torres et al., 2011). The high sequence 
similarity between MERS-CoV and MHVE proteins (Fig. 2C) further 
suggests that the latter also has one a-helical TM domain. 


3.3. MERS-CoV E protein forms pentamers 


We have shown previously that SARS-CoV E protein forms 
homo-pentamers (Li et al., 2014; Parthasarathy et al., 2012). To 
test if MERS-CoVE protein also forms pentamers, Blue-native PAGE 
electrophoresis and sedimentation experiments were performed. 
When MERS-CoV E protein was solubilized in 5mM DPC micelles 
(1:50 protein-to-detergent molar ratio) it migrated as a single band 
in a BN-PAGE gel (Fig. 3A). In common with other membrane pro- 
teins (Gan et al., 2011), the molecular weight of this band may 
not correspond to the molecular weight markers used, or to the 


AQPZ ladder (~25 kDa for the monomer) possibly because of resid- 
ual detergent bound to the protein. However, after dilution with 
25 mM DPC, MERS-CoV E protein dissociated into lower molecu- 
lar weight oligomers, from monomers to pentamers (Fig. 3A). This 
ladder conveniently provides an internal reference that serves as a 
oligomeric size marker. Therefore, by comparison with that ladder, 
we confidently assign the single band observed in 5 mM DPC to pen- 
tameric oligomers. The same oligomeric size was observed when 
using 50 mM DHPC and 20 mM LMPG (Fig. 3A, right panel), which 
are potentially good environments for structural characterization. 

The oligomeric state of E protein was further confirmed by using 
analytical ultracentrifugation. MERS-CoV E protein was solubilized 
in 5mM (C14-betaine micelles at the concentrations described in 
Section 2. Sedimentation equilibrium data could be best fitted to 
a monomer-pentamer self-association reversible model (Fig. 3B 
and C) with apparent association constant, Ka, of 6.3 x 10!7 M~4 
(Log Ka =17.8). This is almost identical to the value (Log Ka =17.5) 
obtained in C14-betaine for the full length wild type SARS-CoV E 
protein in presence of reductant TCEP (Parthasarathy et al., 2012). 
These values are one order of magnitude higher than those obtained 
for truncated Cys-less SARS-CoV E protein (residues 8-65) (Li et al., 
2014) although the latter contained a relatively long (20 residues) 
N-terminal tag, which may have affected the stability of the pen- 
tamer. 

A total of four cysteines are present in the wild-type MERS-CoV 
E protein sequence, two in the predicted TM domain and two in 
the predicted C-terminal juxtamembrane region (Fig. 1). The above 
results suggest that the introduced Cys-to-Ala mutations in the TM 
domain do not affect oligomer stability. Assuming that the mutual 
orientation between the TM helices is the same in all CoV E pro- 
teins, and therefore to that of SARS-CoV E protein (Pervushin et al., 
2009; Torres et al., 2005, 2006), the orientation of the two TM cys- 
teines in MERS-CoV E can be assumed to be near the helix-helix 
interface (Fig. 3D). The TM helix—helix orientation of coronavirus E 
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Fig. 4. Conductance measurements of MERS-CoV E protein in lipid bilayers. (A) 
Histogram corresponding to the current jump amplitudes at +100mV for MERS 
E protein, fitted to a Gaussian model. Experiments were carried out in 500 mM 
KCI solutions at pH 6 in membranes containing 3:1:1 DOPC:DOPE:DOPS. The his- 
togram contains 300 events; (B) MERS E protein current-voltage relationship. This 
experiment is for a single channel, i.e., one experiment which showed almost no 
deviation. 


proteins has been recently validated by in vivo data, where rever- 
tant mutants of SARS-CoV were obtained in mice after introducing a 
channel-inactivating mutation V25F in its E protein. This mutation 
was located at an interhelical position in our pentameric model, 
whereas the mutations that recovered channel activity were at 
positions 37, 30, 26, 27 and 19, all located at the helix—helix inter- 
face and facing V25F (Nieto-Torres et al., 2014). The two cysteines at 
the TM domain of MERS-CoV E are equivalent to positions 30 and 23 
in the SARS-CoV E protein. Despite this putative inter-helical orien- 
tation, these Cys-to-Ala mutations did not seem to affect oligomer 
size or stability. 


3.4, MERS-CoV E protein displays ion channel activity in lipid 
bilayers 


Conductance histograms of recorded traces at +100mV are 
shown in Fig. 4A. We can observe that MERS E protein forms 
an ion channel with conductance of 0.3+0.01nS in 1M KCl. 
The channel displays ohmic behavior (Fig. 4B). These value is 
comparable to those obtained for purified truncated SARS-CoV E 
protein in the same conditions (Li et al., 2014), where conductance 
was 0.39+0.02 nS. Synthetic full length SARS-CoV E protein and 
synthetic transmembrane domain of SARS-CoV E (Ery, residues 
7-38) produced single channel conductances of 0.19 + 0.06 pS and 
0.18 £0.12 nS in 1M NaCl (Verdia-Baguena et al., 2012), although 


these lower values may have been due to modifications introduced 
during synthesis or to impurities. 

In conclusion, we show that MERS-CoV E protein can be 
efficiently purified and, like SARS-CoV E protein, forms homo- 
pentamers where each monomer has a single TM domain. The 
substitution Cys-to-Ala in the two native TM Cys residues seems 
well tolerated and does not affect the pentamer stability or channel 
activity significantly, even though these residues should be located 
at or near the helix—helix interface. The good behavior of this con- 
struct in detergent enables structural characterization that is in 
progress, as well as drug binding biophysical assays. 
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